RM3-TH/01-7 
ROMEl-1326/01 



Next-to-Leading Order QCD Corrections 
to Spectator Effects 
in Lifetimes of Beauty Hadrons 

M. Ciuchini", E. Franco^, V. Lubicz" and F. Mescia^ 



° Dip. di Fisica, Univ. di Roma Tre and INFN, Sezione di Roma III, 
Via della Vasca Navale 84, 1-00146 Rome, Italy 

^ Dip. di Fisica, Univ. di Roma "La Sapienza" and INFN, Sezione di Roma, 
P.le A. Moro 2, 1-00185 Rome, Italy 



Abstract 

Theoretical predictions of beauty hadron lifetimes, based on the heavy quark 
expansion up to and including order l/rn^, do not to reproduce the experimental 
measurements of the lifetime ratios t{B^) /T{Bd) and T(Ab)/r(i?rf). Large correc- 
tions to these predictions come from phase-space enhanced contributions, i.e. 
hard spectator effects. In this paper we calculate the next-to-leading order QCD 
corrections to the Wilson coefficients of the local operators appearing at 0(l/m^). 
We find that these corrections improve the agreement with the experimental data. 
The lifetime ratio of charged to neutral i?-mesons, t{B^)/t{B(i), turns out to be in 
very good agreement with the corresponding measurement, whereas for t{Bs)/t{Bci) 
and T{Ah)/T{Bd) there is a residual difference at the la level. We discuss, however, 
why the theoretical predictions are less accurate in the latter cases. 



1 Introduction 



The calculation of inclusive decay rates of hadrons containing a b quark can be performed 
by computing the amplitude through an Operator Product Expansion (OPE) in powers of 
AQco/^b 01 ■ This allows a quantitative evaluation of the hadronic decay rates under 
the assumption that properly smeared partonic amplitudes can replace the hadronic ones. 
This assumption is usually referred to as the "quark-hadron duality", which in the case 
of total rates is required to hold in its local formulation namely the smearing is 

provided by the sum over exclusive states. 

Within this theoretical framework, and up to terms of 0{l/ml), only the b quark enters 
the short-distance weak decay, while the light quarks in the hadron interact through soft 
gluons only. In particular, the leading term in the expansion reproduces the results of 
the old spectator model, thus providing a theoretical basis to this model. Quantitatively, 
however, the contributions to the lifetimes ratios of the first two terms in the OPE are too 
small. Neglecting terms of 0{l/ml), one finds 

riB+) t(Bs) T(Ab) 

7ky^°-''«- 

to be compared with the experimental measurements @] 

riB+) TiBs) rfAb) 

--— f = 1.07 ±0.02, ^TTT = 0.95 ± 0.04 , ^rrr = 0.80 ± 0.06 . 2 

r{Bd) T{Bd) T{Bd) 

Given the discrepancies, it is important to understand whether we are facing a signal of 
quark-hadron duality violation or, instead, the inclusion of higher-order terms in KQCD/f^b 
and as is sufficient to reproduce the data. 

Spectator contributions, which appear at 0{l/m\) in the OPE, can explain the lifetime 
differences, as they distinguish the light-quark content of the hadrons. As observed by 
Neubert and Sachrajda |^, these effects, although suppressed by an additional power of 
I /nib, are enhanced with respect to leading contributions by a phase-space factor IGvr^, 
being 2^2 processes instead of 1 3 decays. The calculation of spectator effects at the 
leading order (LO) in QCD has been done few years ago [0, and a phenomenological 
analysis has been presented in ref. 0. By using the recent lattice calculations of the 



relevant four-fermion operator matrix elements p-[|12|, we obtain the LO predictions 



t{B^ 



LO 



1.01 ±0.03, 

nBd) 



= 1.00 ±0.01, 

LO 



= 0.93 ±0.04, (3) 

LO 



in agreement with previous estimates. The errors include the uncertainty due to the vari- 
ation of the renormalization scale between mb/2 and 2mf,. For the ratio r(Ab)/r(i?d), 
the inclusion of spectator contributions reduces the discrepancy between theoretical and 
experimental determinations while, for B mesons, at the LO, no significant improvement 
is seen. 

In this paper we compute the next-to-leading order (NLO) QCD corrections to the 
coefficient functions of the valence operators contributing to spectator effects. Our calcu- 
lation has been performed in the limit of vanishing charm quark mass, which means that 
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corrections of the order of asfnl/ml have been neglected. The main motivations to improve 
the leading-order results are 

• reducing the large renormalization-scale dependence which appears at the LO 

• properly taking into account in the matching the renormalization-scheme dependence 
of the Wilson coefficients and of the operator matrix elements; 

• increasing the accuracy of the theoretical predictions by evaluating 0{as) corrections 
which are potentially large. 

Using the results of the NLO calculation, and the available lattice determinations of the 
hadronic matrix elements |0-|12[, we have also performed in this paper a phenomenological 
analysis of the lifetime ratios. We find that the inclusion of NLO QCD corrections to 
spectator effects improves the agreement between theoretical predictions of lifetime ratios 
and their measured values. We obtain the estimates 



t(B- 



1.07 ±0.03. 



NLO 



r{B, 



1.00 ±0.01 



NLO 



t{B, 



0.89 ±0.05, (4) 



NLO 



to be compared with the experimental determinations given in eq. (0). The lifetime ratio of 
charged to neutral 5-mesons turns out to be in very good agreement with the experimental 
data. In the case of the Bs meson and the baryon the agreement is at the la level. 
We should also mention, however, that, in the baryon case, the lattice evaluation of the 
relevant matrix elements is still preliminary . 

An important check of our perturbative calculation is provided by the cancellation of 
the infrared (IR) divergences in the expressions of the coefficient functions, in spite these 
divergences appear in the individual amplitudes. Their presence provides an example of 
violation of the Bloch-Nordsieck theorem in non-abelian gauge theories. We have also 
checked that our results are explicitly gauge invariant and have the correct ultraviolet 
(UV) renormalization-scale dependence as predicted by the known anomalous dimensions 
of the relevant operators. 

The OPE of the lifetime ratios is expressed in terms of local operators defined in the 
Heavy Quark Effective Theory (HQET). The non-perturbative lattice determination of the 
hadronic matrix elements has been performed, in some cases, in terms of operators defined 
in QCD. In order to combine these determinations with the results obtained for the Wilson 
coefficients, consistently at the NLO, a matching between QCD and the HQET operators 
must be performed. In this paper, we have computed this matching, at in the case 

of the fiavour non-singlet AS = four fermion operators, for which the lattice results are 
available so far. 

The NLO predictions of the lifetime ratios may still be improved in several ways. These 
improvements concern both the perturbative and the non-perturbative part of the calcu- 
lation, and are: 

- the charm quark mass corrections, of the order agml/ml, to the Wilson coefficients 
are still unknown. Such a calculation is in progress; 

- the NLO coefficient functions of the current-current operators containing the charm 
quark field and the penguin operator (see eqs. (|l^) and (|l^)) have not been com- 
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puted yet. The non-perturbative determination of the corresponding matrix ele- 
ments is also missing. These contributions do not affect the theoretical determination 
of r(i?+)/r(Srf) and are an S'f/(3)-breaking effect for t{Bs) /r^Bd). In the case of 
T{Ab)/T{Bd), however, their effects may not be negligible; 

the NLO anomalous dimension of the relevant operators in the HQET is still unknown. 



For this reason, present lattice calculations of the Af, matrix elements |T0|, only 
achieve a LO accuracy. For B mesons, a complete NLO determination has been 
performed on the lattice by using operators defined in QCD |T2| ; 



- in the lattice determination of the hadronic matrix elements, penguin contractions (i.e. 
eye diagrams), which contribute in the case of flavour non-singlet operators, have not 
been computed. These contributions cancel in the evaluation of the ratio t{B^) /r^Bd) 
but affect r^Kb) /ri^Bd) and t{Bs) /ri^Bd), the latter through 5'f/(3)-breaking effects. 

The above discussion shows that, while the theoretical determination of t{B^) /r^Bd) can 
be considered as rather accurate, in the case of t{Bs) /r^Bd) and T{Ab)/T{Bd) some ingre- 
dients are still missing, and their determination may further improve the agreement with 
the experimental data. 

The plan of this paper is the following. In sect. ^ we give the relevant formulae for the 
decay rate of 6-hadrons. The details of the calculation are given in sect. |^, while our results 
are summarized in sect. ^. In sect. ^ we present the results of the matching between QCD 
and HQET operators for the flavour non-singlet sector. Finally, the phenomenological 
analysis is presented in sect. |. 



2 The inclusive decay width of beauty hadrons 



Using the optical theorem, the inclusive decay width of a hadron Hf, containing a b quark 
can be written as the imaginary part of the forward matrix element of the transition 
operator T 

T{Hk -^X) = -^lm{H,\r\H,) = ^^{Hb\T\Hb) , (5) 



2M. 



Hb 



where T is given by 



(6) 



The operator ^ is the effective weak hamiltonian which describes AB = 1 transitions. 

It has the form 



n 



AB=1 
eff 



Gp 



V:,Vts E C^Q^ + CsgQsG ] + [s ^ d\ 



+ E + K^QD I + /^-c. , (7) 

l=e,fj,,T ) 
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where Cj are the Wilson coefficients, known at the NLO in perturbation theory |r3|-||15||, 
and the operators Qi are defined as 



Ql = {biCj)v~A{UjSi)v~A, 
Ql = {biCj)v-A{CjSi)v-A, 
Qs = (biSi)v-Aj2(Wj)v-A, 

g 

Q5 = ibiSi)v-A^{qjqj)v+A , 

g 

9s 



Q 



8G 



57r- 



Q2 

Ql 



Ql 



{biCi)v-A{UjSj)v-A , 

{biCi)v-A{CjSj)v-A , 

(biSj)v-A^{qjqi)v-A, 
g 

{biSj)v-A^{qjqi)v+A , 



{biqi)v-A{i^il) 



V-A 



Here and in the following we use the notation {qq)v±A = ql^li^- i Ts)"?- A sum over colour 
indices is always understood. We also neglect Cabibbo-suppressed terms in eq. so that 
the effective Ai? = 1 hamiltonian reduces to 



AS=1 
e// 



\Cj)v-A{Ujd.-^V-A + {l^iCj)v-A{CjSi)v-A 

C2 [{biCi)v-A{Ujdj)v-A + {biCi)v-A{CjSj)v-A 

j2CiQ^ + CsGQ8G+ E Qi]+h.C.. 
i=3 l=e,ii,T ) 



+ 
+ 



(9) 



In the double insertion of the Ai? = 1 effective hamiltonian, this corresponds to neglecting 
terms suppressed by ml /ml sin^ 9c with respect to the dominant ones. 

Because of the large mass of the b quark, it is possible to construct an OPE for the 
transition operator T of eq. (^), which results in a sum of local operators of increasing 
dimension §]. We include in this expansion terms up to 0(l/m^) plus those 1/ml 
corrections that come from spectator effects and are enhanced by the phase space. The 
resulting expression for the local AB = operator F, of eq. (||), is 



192773 



Mlb + c^'^^^ba^,G'"'b 



mt 



mi 



(10) 



where the O^^ are a set of four-fermion dimension-six operators to be specified below. 
These operators represent the contribution of hard spectator effects. At the lowest order 
in QCD, the diagrams entering the calculation of F are shown in fig. |I]. 

We notice that, in general, the OPE of F cannot be expressed in terms of local operators 
defined in QCD. The reason is that, in full QCD, renormalized operators mix with operators 
of lower dimension with coefficients proportional to powers of the 6-quark mass. In this 
way, the dimensional ordering of the OPE is lost. In order to implement the expansion, 
the matrix elements of the local operators should be cut-off at a scale smaller than the 
6-quark mass. This is naturally realized by defining the operators in the HQET. 
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Figure 1: Examples of LO contributions to the transition operator T (left) and to the 
corresponding local operator T (right). The crossed circles represent the insertions of the 
AB = 1 effective hamiltonian. The black squares represent the insertion of a AB = 
operator. 



From eqs. @ and {^^, one can derive an expression for the ratio of inclusive widths 



TjH,) MH'^{bb)H, 



1 + 



(bb) 



+ 



mv 



E 



c(3)(66) 



1 + 



V 



+ 



Hi 



ml 



(11) 



where {■ ■ ■)h denotes the forward matrix element between two hadronic states H defined 
with the covariant normalization. The matrix elements of dimension-three and dimension- 
five operators, appearing in eq. (0), can be expanded using the HQET 



2Mh, 



1 



2ml 



gs{ha^,G^%)H, = 2MH,{2^iUHb) + 0{l/m, 



(12) 
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By substituting these expansions in eq. ([TT|), we finally obtain 




mv, c 



(13) 



which is the expression used in the evaluation of the lifetime ratios of beauty hadrons. 

The Wilson coefficients c*^^-* and c^^-* in eq. ( |T3| ) have been computed at the LO in ref. ||16 
while the NLO corrections to c*^^-* have been evaluated in The NLO corrections 

to c^^-* are still missing. Their numerical contribution to the lifetime ratios, however, is 
expected to be negligible. 

The dimension-six operators in eq. ([T3|), which express the hard spectator contributions, 
are the current-current operators 



with q 



01 = {bb)v-A{qq)v-A , Of = {bb)v+A{qq)v-A , 

01 = {m)v-A{qt''q)v-A , 01 = {bt%)v+A{qt''q)v-A , 
u, d, s, c, and the penguin operator 

Op=ihttjbj)v E (QktMv 

q=u,d,s,c 



(14) 



(15) 



In these definitions, the symbols b and b denote the heavy quark fields in the HQET. 
Note that the operator basis in eq. (|1^ differs from the one used in ref. by a Fierz 
rearrangement. Our choice of the basis is more convenient for NLO calculations because it 
does not require the introduction of Fierz-evanescent operators in the intermediate steps 
of the calculation. 

The coefficient functions of the current operators 0^, have been computed at the LO 



in ref. for q = u,d, s, and in ref. p3| for the charm quark operator. The coefficient 
function of the penguin operator Op vanishes at the LO. 

In this paper we have computed the NLO QCD corrections to the coefficient functions 
of the operators O^ with q = u,d,s. The operators containing the charm quark fields 
contribute, as valence operators, only to the inclusive decay rate of mesons, and their 
contribution to non-charmed hadron decay rates is expected to be negligible. The calcula- 
tion of the NLO corrections to these coefficient functions, as well as the NLO calculation 
of the coefficient function of the penguin operator, has not been performed yet. 



3 Details of the perturbative calculation 

In this section we summarize the general formulae used in the matching procedure, and 
present the details of our calculation. The matching condition between the transition 
operator in eq. (^ and the width operator in eq. (|TD|) can be written in the form 

Im{T)r^c^{f^){d{f,)), (16) 
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where (■ ■ ■) indicates the matrix element computed between a common pair of partonic 
external states. We are using the vector notation for both coefficients and operators. The 
condition (|T6D determines the Wilson coefficients c at the matching scale /i. 

We expand at 0{as) the matrix elements in eq. ([T6|) and express the result in terms of 
tree-level matrix elements (0)o 



T 



QCD 



Im{r) = (To + ^Tij (0)o 
{^{^^)) = {l + ^s){d)^^'\ (17) 
We also consider the perturbative expansion for the Wilson coefficients, 

c = co + ^ci. (18) 

Up to l/irib corrections, the tree-level matrix elements of the operators in QCD are equal 
to their HQET counterparts so that, by using eqs. (|l6])-(0), we readily find the matching 
conditions at the scale /i 

Co = To , ci = fi - To . (19) 



In order to obtain Ti and s, we have first computed in QCD the imaginary part of 
the diagrams shown in fig. |^ {full theory) and then, in the HQET, the diagrams shown in 
fig. ^ {effective theory). The external quark states have been taken on-shell and all quark 
masses, except rrib, have been neglected. More specifically, we have chosen the heavy quark 
momenta pi = ml in QCD and fc^ = in the HQET, and Pq = for all light quarks. In this 
way, we automatically retain the leading term in the l/mh expansion. We have performed 
the calculation in a generic covariant gauge, in order to check the gauge independence of 
the final results. Two-loop integrals have been reduced to a set of known (massless p- 
and massive tadpole) integrals using the recurrence relation technique Equations 
of motion have been used to reduce the number of independent operators. The AB = 
operators in the HQET have been renormalized in the NDR-MS scheme defined in details 



in ref. 27 



Some diagrams, both in the full and in the effective theory, are plagued by IR diver- 
gences. These are diagrams where a soft gluon is exchanged between the external legs {Di, 
D2, D^, Dq in fig. ^ and all the diagrams in fig. i Q. The diver gences do not cancel in 
the final partonic amplitudes, and provide an example of violation of the Bloch-Nordsieck 
theorem in non-abelian gauge theories [Hf-I^. In agreement with this theorem, we explic- 
itly verified that the abelian combination of the diagrams does not contain IR divergences. 
They only appear when the colour structure is taken into account. According to the KLN 



theorem [^, |3^, IR singularities cancel when the contribution of soft gluons in the initial 
state is included in the amplitudes. Even if the amplitudes in the full and effective theo- 
ries are IR divergent, we expect the IR poles to cancel in the matching, since the Wilson 
coefficients must be insensitive to soft physics. We have explicitly checked that, in the 
computation of the coefficient functions c^^\ this cancellation takes place. 

^ We notice that the introduction of non-vanishing masses and momenta for the external Hght quarks 
is not sufficient to eUminate the IR divergences. 
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Figure 2: Feynman diagrams which contribute at NLO to the matrix element of the tran- 
sition operator T in the case q — s. In the other cases, q — u, d, diagrams D13 and D14 
do not contribute. 
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Figure 3: Feynman diagrams which contribute at NLO to the matrix element of the AB = 
operators entering the width operator T. 

We use D-dimensional regularization with anticommuting 75 (NDR) to regularize both 
UV and IR divergences. The presence of dimensionally-regularized IR divergences intro- 
duces subtleties in the matching procedure. Usually, the matching can be performed by 
only considering the four dimensional operator basis, since renormalized evanescent oper- 
ators do not give contributions to the physical amplitudes. In the present case, however, 
IR poles in e = (4 — D) /2 promote the 0{e) contribution from evanescent operators to 
finite terms in the D ^ 4 limit. This contribution has to be taken into account and the 
matching procedure must be consistently performed at 0{e). To illustrate this point in 
further details, we provide the example in fig. ^. Let us consider the IR-divergent parts 
of the diagrams Di and Ei. These divergences are not expected to give contribution to 
the Wilson coefficients, since the coefficients should only account for the UV behaviour of 
the effective theory. This condition is in fact guaranteed by the factorization of IR diver- 
gences and by the LO matching condition. Indeed, the numerator entering the definition 
of A/ijCf^Lo fig- @5 which represents the IR contribution to the coefficient function c^^-* 
at the NLO, vanishes just because of the LO matching condition. This happens provided 
that IR divergences are regularized in the same way in the full and in the effective theory, 
independently of the specific choice of the IR regulator. In this way, soft physics does not 
enter the coefficient functions. If, however, IR divergences are dimensionally regularized, 
the condition A/ijcj^Lo ~ only holds when the LO matching condition is performed up 
to and including terms of 0{e). In particular, D-dimensional matching requires enlarging 
the operator basis to include (renormalized) evanescent operators, which must be inserted 
in the one-loop diagrams of the effective theory. Because of the IR divergences, the matrix 
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Figure 4: Example of the mechanism that guarantees the vanishing of IR contributions to 
the coefficient functions. Only the contributions from diagrams Di and Ei are shown. In 
the presence of dimensionally-regularized IR divergences, LO matching up to and including 
terms of 0{e) is required. 



elements of the renormalized evanescent operators do not vanish in the D — > 4 hmit. 

A similar discussion concerning the matching procedure in the presence of IR diver- 
gences can be found in ref. ||3^. We mention that the relevant evanescent operators, which 



enter the matching in our calculation, are 

Ef = ht^lalu Lbq^TlU-^Ol, El = hf.la'lu it^b qYTllf'q -^Ol, 

(20) 

El = h^-i^-i,Rb qrrilq - 16 01 , El = h^laluRf'b qYl'^llt^q - 16 0| . 

The procedure for taking into account their contribution is further discussed in the Ap- 
pendix. We have also checked, on a proper subset of gauge-invariant diagrams, that a 
calculation using the gluon mass as IR regulator gives the same result as the one obtained 
with dimensional regularization by performing the matching at 0{e)^ and including the 
insertion of the evanescent operators. 

We conclude this section by mentioning that our calculation of the Wilson coefficients 
has passed the following checks: 
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• gauge invariance: we verified that the coefficients are exphcitly gauge- invariant. The 
same is true for the full and the effective amplitudes separately; 

• renormalization-scale dependence: the coefficient functions have the correct renorma- 
lization-group behaviour as predicted by the LO anomalous dimension matrix of the 
AB = operators; 

• IR divergences: the coefficient functions are infrared finite. The cancellation of IR 
divergences also takes place for the abelian combination of diagrams in both the full 
and the effective amplitudes. 



4 Spectator contributions at the NLO 

In this section we present the results for the coefficient functions of the dimension-six 



operators Ol of eq. ( [T^ ) for q = u,d, s. We collect the known LO expressions and then 
present our new results for the NLO contributions. 

The width operator F in eq. (|^) depends quadratically on the coefficient functions Cj 
of the AB = 1 effective hamiltonian. Therefore, we find convenient to write the coefficient 
functions of the dimension-six operators as 

4(/^)= E (21) 

where fii is the renormalization scale of the AB = 1 effective hamiltonian and fi is the 
renormalization scale of the AB = operators entering the expansion of F 0. The coef- 
ficients Fj}- depend on the renormalization scheme and scale of both the AB = and 
AB = 1 operators. The dependence on ni and on the renormalization scheme of the 
AB = 1 operators actually cancels, order by order in perturbation theory, against the 
corresponding dependence of the AB = 1 Wilson coefficients Cj. Therefore, the coefficient 
functions only depend on the renormalization scheme of the AB = operators. A given 
scheme of the AB = 1 operators, however, must be chosen in order to present results for 
F^-j. In the following we always consider AB = 1 operators renormalized in the NDR-MS 
scheme, defined in details in ref. MM. The corresponding Wilson coefficients can be found 



in refs. 13 - 15 



Since we use dimensional regularization to regularize both UV and IR divergences, the 
scales fi and /ii in eq. (^) play the role of both renormalization scale and IR regulator. 
Dimensionally-regularized IR divergences produce additional log/x and log/xi, beside those 
governed by the renormalization-group equations. In the matching procedure, the same IR 
regulators in the full and effective theories must be used, therefore we have to identify fi 
and /ii. This choice gets rid of the spurious log(/i/yUi) so that the NLO coefficient functions 
have the correct UV behaviour predicted by the LO anomalous dimensions of the AB = 
operators. 

For the coefficients F^ -j we write the expansion 
/^) = + ^BIM) ■ (22) 



^ In the case q — s, at the NLO, there are other contributions to eq. ( pi| ) coming from the insertion of 
AB = 1 penguin operators > 2) which wiU be discussed at the end of the section. 



11 



Since, by definition, the coefficients F^-- are symmetric in the indices i and j, we will only 



present results for i < j. The leading order coefficients Af- - have been computed with 



k,ij 

a non- vanishing charm quark mass by Neubert and Sachrajda [j^. The relevant one-loop 
diagram is shown in fig. ^ (left). We have repeated the calculation and found results in 
agreement with ref. [0. In the basis of eq. ([l^), the coefficients Al ^j, for the case q = d, 
read 



3 V 2/ ' ' ' 6 
^li2 = -^(l-^)' (1 + 2^) , = -^(1-^)' (1 + 2^) , 

4ii = -2 (l + (1 - , 4i2 = (l + i) (1 - ' 



(23) 



^22 = 0, = -(1-2;)' (1 + 2 2:) , 

1 

3' 



<12 = -5(1 -^)' (1 + 2 2) , <22 = 0. 



where z = m^/m^. 

At the LO, the only difference between the case q = s and q = d is the massive charm 
quark running in the loop. The coefficients Af. are given by 

Al,, = -^VT^ {1-z), Al,, = -^VT^ {l~z), 

1 , , . .„ 1 



^22 = -3^1^(1-^), Al,, = --Vl^z {l + 2z) 

1 , , . I 



4i2 = -T^Vl^4^ (1 + 2;^) , = -7Vr^4i {l + 2z) 

18 ' D 



(24) 



Al,, = -2Vl^z {l-z), Al,^ = --VT^z {l-z), 

K22 = , 4ii = -v^l^ (1 + 2 z) , 

A|_,2 = (1 + 2 z) , A|^22 = . 

For q = u instead, the relevant one- loop diagram is shown in fig. H (right). It gives 



^1,11 


= (l-^)% 


AU 

^1,12 




4" 

^1,22 


= {l-zf 


AU 

^2,11 


= 0, 


AU 

^2,12 


= 0, 


AU 

^2,22 


= 0, 


AU 

^3,11 


= 0, 


AU 

^3,12 


= 2{l-z)\ 


AU 

^3,22 


= 0, 


AU 

^4,11 


= 0, 


AU 

^4,12 


= 0, 


^4,22 


= 0. 



(25) 
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Figure 5: Feynman diagrams entering at the LO in the calculation of the coefficient func- 
tions cl for q = u,d,s. 



The NLO coefficients B'^-j are obtained by computing the diagrams shown in figs. |^ 
and 0. The details of the matching procedure have been discussed in the previous sec- 
tion. We just remind here that we have neglected corrections of 0{ml/ml) in the NLO 
calculation. 

At the NLO, the coefficients depend on the renormalization scheme chosen for the 
HQET operators in eq. ([T^) . We have computed these coefficients in the NDR-MS scheme 
defined in ref. [^. In this scheme, for the coefficients B^^j{fi) we obtain 

_ 230 _ 778 32L _ 314 , 16n' 



81 ' '''^ 243 ' 9 ' ''^^ 81 27 ' 

_ 404 16L _106 Svr^ 

^2,11 - , ^2,12 - 243 + 9 ' ^'22 " 81 ^27 ' 

-Bq 1 1 = h 1 8 , i?o 1 9 = 1 , B'i 99 = h 8 

3,11 54 ' 3,12 3 ^'^^ 9 9^ 

Ba 1 1 = h 9 , B'i 1 9 = 1 , B'i 99 = h 4 Zy 

4,11 54 ' 4,12 3 ^'^^ 18 9 



(26) 



where L = log {fi/rrib). In the case q = u, the relevant Feynman diagrams involve gluon 
corrections to the LO diagram shown on the right side of fig. S. We obtain for the coefficients 



TDU 

-Dl,ll - 


188 


16 TT^ 

9 ' 




TDU 

-"1,12 - 


1312 

81 


32 7r2 
27 


r)u _ 

-°1,22 ~ 


188 


16 7r2 
9 ' 




TDU _ 

-°2,11 — 


16 
27' 




TDU 

-"2,12 ~ 


32 
~81 ' 






Oil 

-°2,22 ~ 


16 
27' 




TDU _ 

-°3,11 ~ 


125 


4 7r2 


-24L, 


TDU 

-'-'3,12 ~ 


374 


28 7r2 


3 


3 


27 


9 


TDU 

-^3,22 ~ 


125 
3 


4 7r2 
3 


- 24L, 


DM 

-D4,ll - 


4 
~3 ' 




TDU _ 

-°4,12 — 


100 
27 






TDU 

-^4,22 ~ 


4 
~3 ■ 





32 L 
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g = 


= d 


1 = 


= u 


q = 


= s 




LO 


NLO 


LO 


NLO 


LO 


NLO 


cl 


-0.32 


-0.30 


0.92 


0.90 


-0.27 


-0.25 


cl 


-0.18 


-0.16 


0.00 


0.02 


-0.17 


-0.15 




0.20 


0.16 


-0.83 


-1.42 


0.17 


0.13 


C4 


0.11 


0.08 


0.00 


0.00 


0.11 


0.08 


Cl 


-0.02 


-0.03 


-0.06 


-0.33 


-0.01 


-0.03 


C2 


0.02 


0.03 


0.00 


-0.01 


0.02 


0.03 


Cz 


-0.70 


-0.65 


2.11 


2.27 


-0.60 


-0.54 


C4 


0.79 


0.68 


0.00 


-0.06 


0.77 


0.65 



Table 1: Wilson coefficients cl computed at the LO and NLO (MS scheme). As reference 
values, we use /i = rrih = 4.8 GeV and rric = 1.4 GeV. The coefficients cj, in the operator 
basis of eq. are also shown. 



Finally, we discuss the extension of eq. (^) to the case q = s. With respect to cf.i^fi), 
the coefficient functions c|(/i) receive, at the NLO, additional contributions coming from 
the penguin contraction of the current-current operators (diagram D13 of fig. |^) and from 
the insertion of penguin and chromomagnetic operators (diagram D14 of fig. |^). Since 
the Wilson coefficients C^-Cq are small, contributions with a double insertion of penguin 



operators can be safely neglected. As suggested in |^5[, a consistent way for implementing 
this approximation is to consider the coefficients C^-Cq as formally of 0{as)- Within this 
approximation, only single insertions of penguin operators need to be considered at the 
NLO. Therefore we can write 



j=l,2r=3,4 



(2^ 



and we obtain 



-Pi, 22 

-P3,22 



32 
243 



32 L 



81 
? 8L 

1 ~ 27 



-P2,22 
-P4,22 



+ 



16 L 



16 

243 ' 81 
4 4L 

~81 ~ '27 



(29) 
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-Pi, 13 = 


1 

~3' 


-P2,13 = 


1 

^6' 


-P3,13 = 


-2, 


-P4,13 


-Pi, 14 = 


1 

~9' 


-P2,14 = 


1 

"18 ' 


-P3,14 = 


2 

"3' 


-P4,14 


-Pi, 23 = 


1 

~9' 


-P2,23 = 


1 

~18' 


-P3,23 = 


2 

~3' 


-P4,23 


-Pi, 24 = 


1 

~3' 


-P2,24 = 


1 

~6 ' 


-P3,24 = 


0, 


-P4,24 


-Pi, 28 = 


8 

"27' 


-P2,28 = 


4 

~27' 


-P3,28 = 


2 

9' 


-P4,28 



-1 



1 

■3' 

1 

■3' 



(30) 



0, 
1 

9 ■ 

The coefficients Pfc,28 are computed using the convention in which the chromomagnetic 
coefficient Cg has a positive sign. The NLO contribution of penguin and chromomagnetic 



operators to beauty hadron hfetimes has been also computed in ref. ||3^. We verified that 
our results are in agreement with them. 

For convenience, we present in table |1| the numerical values of the coefficients c^(/i), 
in the cases q = u,d, s, at fi = nib = 4.8 GeV, both at LO and NLO. We also give the 
coefficients c^(/u) defined in the operator basis of ref. 0]. For details on the conversion 
between the two basis, see sect. BTTI 



5 Matching of AB = non-singlet operators between 
QCD and HQET 

Theoretical predictions of beauty hadron lifetimes, based on the OPE, are obtained by com- 
bining the perturbative determination of the Wilson coefficients with the non-perturbative 
evaluation of the hadronic matrix elements of the HQET operators. In the next section, 
we will perform a phenomenological analysis, by using the NLO results for the Wilson 
coefficients and the lattice determinations of the relevant hadronic matrix elements. In 
some cases the lattice results are provided in terms of matrix elements of QCD operators. 
In order to use these determinations, consistently at the NLO, it is necessary to perform 
an additional matching between the QCD and the HQET operators. For the dimension-six 
AB = operators of interest in this paper, this matching has not been computed so far. 
For this reason, we present in this section the results of an 0{as) calculation of the coeffi- 
cient functions relating the QCD four-fermion AB = operators to their counterparts in 
HQET. 

We limit ourselves to the case of fiavour non-singlet operators, which are the only 
ones for which lattice results are available so far. An important example is the difference 

— of of the operators defined in eq. (0) which determines, as we will see below, the 
lifetime ratio of charged to neutral B mesons. In the limit of exact fiavour symmetry, these 
operators do not mix with lower dimensional operators. For this reason, their matching 
between QCD and HQET does not require the calculation of penguin contractions and 
only involves current-current operators of dimension six. 

The matching equation between the operators in QCD and HQET can be written in 
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the form 

OQCD(/i) = <5(Ai, /i', rrih) OnqETifJ'') ■ (31) 

The scale dependence of the coefficient functions C{fi, fi' ,mb) is governed by the renormal- 
ization group equation of the QCD and HQET operators respectively. Thus one finds 



/i', rrih) = {W{n, rribf^ ^ C{mi,) W{^i\ m^f 
where the matrix W^(yU, mb) in QCD is given, at the NLO, by 



(32) 



1 + '^J{^^)\ { { 1 - ^ j(m,) 



47r 



An 



(33) 



and W{n,mii) is its analogous in the HQET. For convenience, we present here the expres- 
sions of the leading anomalous dimension matrices, 70 and 70 in QCD and HQET. In the 
basis of eq. (|ll they reads IHEI 



7o 



/ 








12 





\ 













-12 






8/3 





-4 







V 





-8/3 





-14 


/ 



7o 



( ° 

































-9 







lo 








-9 


/ 



(34) 



The matrix J, in eq. (p3D, is expressed in terms of the two-loop anomalous dimension, 
and it is defined for instance in ref. p4|. In the case of the HQET, the NLO anomalous 
dimension is still unknown. 

The matrix C(mf,) is the coefficient function at the matching scales /i = /i' = m^. At 
the NLO, it can be written in the form: 



C[nib) = IH Ci . 

47r 



(35) 



We have computed Ci considering two different renormalization schemes for the QCD 



operators, namely NDR-MS and Landau RI-MOM (see ref. for a detailed definitions of 
these schemes). The HQET operators are renormalized in the NDR-MS scheme of ref. [2? . 
By denoting the results with C^'^ and C^^ respectively, in the operator basis of eq. ([l^), 
we obtain: 

/ 



MS 



-4/3 


4/3 


-18 


-2 \ 


4/3 


-4/3 


2 


6 


-4 


-4/9 


-35/6 


-5/2 


4/9 


4/3 


-5/6 


25/6 y 



(36) 



■^Note that the expression of the LO anomalous dimension matrix in HQET given in eq.(A.l) of ref.Q] 
contains a misprint: the matrix element (70)12 should be read as -1 instead of 1. 
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and 



-4/3 
4/3 

V9- 16/3 log 2 
4/9 



4/9 



4/3 
-4/3 
-4/9 

-8/9 log 2 



-4 - 24 log 2 
2 

-21/2 + 8 log 2 
-5/6 



-2 \ 

2-41og2 

-5/2 
-13/2 + 4/3 log 2 J 



(37) 



Note that the difference Cf^^ — Cf-^ is the one loop matrix which provides the connection 
between the MS and RI schemes in QCD for the AB = operators. We have checked that 
this matrix agrees with the matrix denoted as in ref. |^|, once the proper change of 
basis is performed. 

In the next section, the matrix Cf^^ of eq. (|36| ) will be used to convert the lattice results 
for the matrix elements of the QCD operators into their HQET counterparts. 



6 Phenomenological Discussion 

In this section we perform a phenomenological analysis of the lifetime ratios of beauty 
hadrons, by using the NLO expressions of the Wilson coefficients and the lattice determi- 



nations of the relevant hadronic matrix elements 0-112]. 

The starting point of this analysis is eq. (|I3]), which expresses the ratio of inclusive 
widths of beauty hadrons, up to and including 1/ml spectator effects. The combinations 
of hadronic parameters entering this formula at order I /ml can be evaluated from the 



heavy hadron spectroscopy |39]. The numerical contributions of these terms to the lifetime 



ratios are rather small, and one obtains the estimate 

^^ = 1.00-Afp:,, li|£l = l.OO-Afp;,, lM = 0.98(1) -Aj,,,, (38) 

T{I^d) 'T{^d) '^{^dj 

which can be compared with the experimental results in eq. (^ . In the previous expressions 
the As represent the 1/ml contributions of hard spectator effects 

These are the quantities which we are interested in. They are expressed in terms of coeffi- 
cient functions and matrix elements of dimension-six operators. In the phenomenological 
analysis, we use the non-perturbative determination of the relevant matrix elements from 



lattice QCD calculation. For recent estimates based on QCD sum rules, see refs. pO[-|^|. 
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Figure 6: Feynman diagrams entering the non-perturbative evaluation of the matrix ele- 
ments of valence operators between external B-meson states. In the case of non-valence 
operators only the second diagram appears. 

6.1 Matrix Elements 

The basis of four-fermion operators usually considered in the literature to study the lifetime 
ratios of beauty hadrons is given by [0 



Ol = {hiqi)v~A {qjbj 



'jUj)V-A , 



Ol = {ht1^q,)v-A (qkttik 



'l)V-A 



Ol = ipiqi)s-p {qjbj)s+p, 

Ol = {ht%q,)s^P {qktMs+P, 



(40) 



where q = u,d,s,c and {qq)s±p = q{^ =t 7^)13'. In this basis, which differs from the one 
in eq. by a Fierz rearrangement, the matrix elements of the operators and Of, 
between external 5-mesons states, vanish in the vacuum saturation approximation (VSA). 
By following the convention adopted in the literature, we will work in this section in the 
basis of eq. pO|) , to which we add the penguin operator of eq. (|15D . The Wilson coefficients 
associated to the new basis will be indicated by the symbol c^. They are related to those 
defined in the previous sections by a linear transformation 



M 



1 


1/3 





4/9 












-2/3 





-8/9 






2 





-1/3 












-4 





2/3 


/ 



(41) 



The values of at /i = m;, are given in table I. 

In parametrizing the matrix elements of the current-current operators, it is useful to 
distinguish two cases, depending on whether or not the hght quark q of the operator 
enters as a valence quark in the external hadronic state. From a diagrammatic point of 
view, this difference gives rise to a different number of Wick contractions. In the case of 
external i?-meson states, for instance, the matrix elements of the valence operators are 
computed by evaluating non-perturbatively the two Feynman diagrams shown in fig. 
while for non-valence operators only the second diagram appears. We find it convenient to 
introduce different i?-parameters for the valence and non-valence contributions. Thus, for 
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the 5-meson matrix elements of the non-valence operator we define 



■^Sr forg^g', (42) 



2Mb 2 
and, for the valence operators (g = g'), we write 



is'? + 6?''] 



'^Mb, 2 ^ ' ' 2Mb, 2 

{B^O^_flMB^ (AM\Bl_flMB^ 



In eq. (|43|), the d^'^ are defined as the parameters 5^'^ of eq. ( ^2]) in the limit of degenerate 
quark masses (m^ = ntqi). In the VSA, Bi = B2 = 1 while the e parameters and all the 
6s vanish. It should be clear that the B^, Sk and 6k have a direct interpretation in terms of 
Feynman diagrams, and are represented by the first {Bk, Sk) and second {6k) contraction of 
fig. ^ respectively. Each of these parameters is a well defined and gauge invariant quantity. 
In the exact SU{2) limit, for instance, the parameters Bj^ express the matrix elements of 
the non-singlet operator — between external 5-meson states. 

Our definition of Bi, B2, Si and 62 differs from the one considered in the literature for 
the presence of the 6s in eq. The reason why we have distinguished between valence 



and non valence contributions, in these definitions, is that lattice QCD calculations P|-||12 
have not computed, so far, the full matrix elements, but only the valence contributions 
represented the first diagram of fig. ^. Thus, the Bk and Sk parameters, defined as in 
eq. (|43|) , denote the quantities actually determined, at present, by lattice calculations. In 
principle, a non-perturbative calculation of the 6 parameters on the lattice is possible. 
However, it requires dealing with the difficult problem of power-divergence subtractions 
which has prevented so far the calculation of the corresponding diagrams. 

To complete the definitions of the i?-parameters for the i?-mesons, we introduce a 
parameter for the matrix element of the penguin operator 

{B,\Op\Bl^flMBp,_ ^^^^ 



2M 



Ba 



We now define the 5-parameters for the A;, baryon. Up to 1/mf, corrections, the matrix 
elements of the operators O2 and 0\, between external A;, states, can be related to the 
matrix elements of the operators Of and Ol [0 

{K,\Ol\k,) = -2 {A,\Ol\k,) , {MOl\K) = -2 (A,|0||A,) . (45) 

For the independent matrix elements, assuming SU{2) symmetry, we define 

2Ma, 2 V ^ ^ ^ 

(A^I^IIA.) PbMb ^^A,x , , 
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Bf 


= 1.2 ± 0.2 


Bf-- 


= 1.0 ± 0.1 


Bl 


= 0.7 ± 0.1 




= 0.7± 0.1 


el = 


= 0.03 ± 0.02 


= 


0.03 ± 0.01 


4 = 


= 0.04 ± 0.01 


^2 = 


0.03 ± 0.01 


Li - 


= -0.2 ±0.1 


L2 -- 


= 0.2 ±0.1 


rrib = 


4.8 ±0.1 GeV 


mb-rric = 


3.40 ± 0.06 GeV 


fB = 


200 ± 25 MeV 


fss / fs 


= 1.16 ±0.04 



Table 2: Central values and standard deviations of the input parameters used in the nu- 
merical analysis. 



2Ma 



(46) 



2Ma, 
{K\Op\^b) 



— O2 tor g = s, c . 



2Ma 



In analogy with the 5-meson case, the parameters Li and L2 represent the valence contri- 
butions computed by current lattice calculations |T^, |TT|. 

At present, two independent lattice calculations of the ^-parameters in the S- meson 
sector have been performed, both in the quenched approximation. In the first study P, p 



the parameters -Bf, -Bf, and 82 have been computed by simulating on the lattice the 



HQET. The second calculation [|12|, instead, has been performed by using QCD, and the 
results have been obtained in this case for both Bud and Bg mesons. The calculation 



of ref. |]12| is also accurate at the NLO, since the -B-parameters have been evolved, from 
the lattice scale up to mf,, by using the two-loop anomalous dimension of the four-fermion 
A-B = operators, which is known in QCD 0, but not in the HQET. For these reasons, 
we will use the results of ref. [12| in our phenomeno logical analysis. The parameters -B^, 



-B2, £1 and computed in QCD, have been matched from QCD to HQET, at the NLO, 



using the coefficients presented in the previous section, eqs. (|35D -(p6D. After applying these 
coefficients to the QCD lattice results of ref. [|12|, we obtain the estimates of the HQET 
-B-parameters collected in table ||. 

For comparison, we also present here the values of -B-parameters obtained in ref. [11] 
from the lattice simulation in the HQET. In this different definition of the MS 

scheme has been adopted. The NLO connection between this scheme and the one of 
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ref. chosen in this paper, is given by [^] 



Bl 
el 



1 + ^ 
47r I -1/3 




Bl 

4 



(47) 

SDP 



where the label SDP denotes the parameters in the scheme of ref. 0. The other param- 
eters, Bl and eX, are the same in both schemes. After applying eq. (^7|), the results of 



refs. P, [Tl| read 



5f = 1.06 ±0.08, 5^ = 1.01 ±0.07, 

e^f = -0.01 ±0.03, = -0.03 ±0.02. (48) 

The differences between these results and those given in table ^ may give an estimate of 
the higher orders l/m^, corrections neglected in the static calculation of ref. |^, |ll|], and 
of the systematic uncertainty introduced by the large extrapolation in the 6-quark mass 
performed in ref. It should be also kept in mind, however, that the renormalization 



group evolution has been only performed in ref. [Tl| with a LO accuracy. 

For the baryon, only the non-perturbative results of an exploratory study in the 
HQET are available at present [0. They have been obtained with a LO accuracy, at 
a rather large value of the lattice spacing and do not include the extrapolation of the 
light quark masses to their physical values. For these reasons, in quoting the values of the 
corresponding B parameters, Li and L2 in table ^, we also include in the error our estimate 
of the remaining systematic uncertainties. We also note that, in all these calculations, the 
matching between the lattice and the continuum theory has been performed using 1-loop 
perturbation theory, and that the leading perturbative corrections are typically found to 
be large. For this reason, a non-perturbative evaluation of the relevant renormalization 



constants would be very interesting. For B mesons, such a calculation is in progress . 

For the large number of 5 and P parameters in eqs. (^2D-(^UD, which define the non- 
valence contributions and the matrix elements of the penguin operator, only phenomeno- 
logical estimates exist P5| , ^7\ . These estimates indicate that the matrix elements of non- 
valence operators are suppressed, with respect to the valence ones, by at least one order 
of magnitude. These contributions cancel out in the theoretical evaluation of the lifetime 
ratio T{B~^)/T{Bd), see eq. (^0]), while they enter the determination of r(A;,)/r(i?rf) and 
t{Bs) /T{Bd)i the latter for 5't/(3) breaking effects. Lacking quantitative calculations of 
the 5 and P parameters, we will mainly rely in our numerical analysis on the results of 
the phenomenological estimates, and neglect these contributions. A qualitative estimate 
of the error introduced by this approximation in the case of the ratio T{Ab)/T{Bd) will be 
given at the end of the section. 



6.2 Results 

We now present the detailed expressions of the spectator contributions, represented, in 
eq. (^9|), by the quantities Agpec- In these expressions all terms, coming from both valence 
and non-valence contributions, which are often omitted in the literature, will be explicitly 
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taken into account. Using eq. (P^D, and the definitions of tlie 5-parameters in eqs. (J^)- 
(^), tfie spectator contributions take tlie following expressions 



AB+ aq^2 Ib^B (~u ~d\ Kid 



"^feC^ U=i 

5r - St) + - 5.^') ] + CP (r P ^ - P } , (49) 

(4^^ - St) + ~ct [dt - St)] +cp[P^-P')}. 



The Cfc are the coefficient functions given in eq. The factor r denotes the ratio 

{fss^Bs) / Ub^b) and, in order to simphfy the notation, we have defined the vectors of 
parameters 

& = {BlBlele\}, 

C = {L,,-Lj2,L2,-L2/2}, (50) 

Because of the SU{2) symmetry, the non- valence and penguin contributions cancel 
out in the expressions of the lifetime ratio t{B^)/t{B(]). From this point of view, the 
theoretical prediction of this ratio is at present the most accurate, since it depends only on 
the non-perturbative parameters actually computed by current lattice calculations. The 
prediction of the ratio r(Af,)/r(Prf), instead, is affected by both the uncertainties on the 
values of the 6 and P parameters, and by the unknown expressions of the Wilson coefficients 
Cfc and cp at the NLO. For the ratio t{Bs) /r^Bd) the same uncertainties exists. However 
they are expected to be smaller, since these contributions cancel in the limit of exact SU (3) 
symmetry. 

In order to evaluate the lifetime ratios, we have performed a Monte Carlo calculation, 
by extracting the input parameters with flat distributions with central values and standard 
deviations given in table |^. The contributions of all the 5 and P parameters have been 
neglected. The strong coupling constant has been kept flxed at the value ^^(mfe) = 0.214. 
The values of the bottom and charm quark masses, given in the table and used in the 
numerical analysis, correspond to the pole mass deflnition. The parameter c*^^^ in eq. ( pU] ) 
is a function of the ratio ml/ml, and such a dependence has been consistently taken into 
account in the numerical analysis. For the range of masses given in table |^, c^^^ varies in 
the interval c^^) = 3.4 ^ 4.2 |2D|-|2g. 



For convenience, we also present the numerical expression corresponding to eq. (pOD, by 
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omitting the non valence and penguin contributions. We find 



spec 



spec 



spec 



- 0.011(3) Bf - 0.005(1) + 0.7(2) ef - 0.19(5) , 



- 0.004(1) B'^ + 0.005(1) Bl - 0.17(4) + 0.22(6) £^ 
+ 0.004(1) Bf - 0.005(1) Bi + 0.17(4) ef - 0.19(5) ef, ^ LO 



- 0.022(6) Li + 0.24(6) La 

+ 0.004(1) Bf - 0.005(1) + 0.17(4) ef - 0.19(5) ef , 



(51) 



at the LO {/i = nib), and 



spec 



spec 



A^ 

spec 



0.07(2) Bf - 0.011(3) Bf + 0.7(2) - 0.18(5) e 



d 
2 ) 



- 0.007(2) Bl + 0.009(2) 5| - 0.16(4) + 0.20(5) 4 
+ 0.007(2) Bf - 0.008(2) 5^ + 0.16(4) ef - 0.16(4) £^ , 

- 0.09(2) Li + 0.32(8) L2 

+ 0.007(2) Bf - 0.008(2) 5^ + 0.16(4) ef - 0.16(4) £^ , 



\ NLO 



(52) 



at the NLO. Note that we have shown, for each coefficient, an estimate of the error. These 
errors, however, are strictly correlated. For this reason, eqs. (|5T| ) and (0) have not been 
used in the numerical analysis. 

We finally present our predictions for the lifetimes ratios of beauty hadrons, which have 
also been quoted in the introduction. Using the LO expressions of the Wilson coefficients, 
we obtain 



r{B, 



1.01 ±0.03, 



LO 



r{Bs 
r{Ba 



1.00 ±0.01 



LO 



r(B, 



0.93 ±0.04, 



(53) 



LO 



while, including the NLO corrections computed in this paper, we get 



r(5^ 



1.07 ±0.03 



NLO 



r{B, 
r{B, 



1.00 ±0.01 



NLO 



r(Afe 
riB, 



0.89 ±0.05. (54) 



NLO 



The central values and errors quoted in eqs. (^3|) and (^4|) are the average and the standard 
deviation of the theoretical distributions. These distributions are shown in fig. |^, together 
with the experimental ones. 

In the LO results, the errors include the effect of varying the renomalization scale fi 
between rub/^ and 2mb. At the NLO, the variation with fi of the Wilson coefficients is com- 
pensated by the fi dependence of the renormalized operators, up to NNLO corrections. We 
cannot estimate the residual scale dependence because we only know the non-perturbative 
the relevant operators at a single scale /i = rrib. Therefore the renormalization scale is kept 
fixed in the NLO analysis. 
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Figure 7: Theoretical (histogram) vs experimental (solid line) distributions of lifetime ra- 
tios. The theoretical predictions are shown at the LO (left) and NLO (right). 
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We find that, with the inclusion of the NLO corrections, the theoretical prediction 
for the ratio T{B'^)/T{Bd) turns out to be in very good agreement with the experimental 
measurement, given in eq. (^. For the ratios t{Bs) /r^Bd) and T{Ab)/T{Bd) the agree- 
ment is also very satisfactory, and the difference between theoretical and experimental 
determinations is at the la level. 



The errors quoted in eq. (5^) do not take into account the systematic uncertainty due to 
terms of 0{as'm'^/'ml) ~ 0.1 a^, not included in our calculation. Naively, these corrections 
can change the NLO terms by about 10%, thus affecting the lifetime ratios at the level of 
1%. A realistic estimate, however, requires an explicit calculation, which is under way. 

We conclude this analysis with a qualitative estimate of the contribution of the non 
valence and penguin matrix elements, which have been neglected in the evaluation of 
t{Bs) /r^Bd) and T{Ah)/T{Bd). We consider as an example the term 



2 4 



48vr^^E[(c7 + c.^)K^-^f)] (55) 



entering the expression of A^^^, in eq. (BO). Since the coefficient functions and are 



known at the NLO, the only uncertainty in eq. ( ^5]) is the value of the non-perturbative 



parameters. According to the phenomenological estimates of ref. I^^l^ typical values 



of these parameters are smaller by approximately one order of magnitude with respect to 
the matrix elements of the valence operators. By choosing for all the differences c?^'' — ^^"^ 
a common value of 0.05, we then find that the ratio r(Af,)/r(i?rf) varies by less than 3%, 
and the variation is approximately proportional to the values of the S parameters. The 
matrix elements of the penguin operators are not expected to be smaller than those of 
the valence operators. Numerically, since the coefficient function cp vanishes at the LO, 
this contribution is expected to have the size of a typical NLO corrections. It is clear 
that a quantitative evaluation of the non-valence and penguin operators would be very 
interesting, in order to further improve the agreement between theoretical and experimental 
determinations of the lifetime ratios of beauty hadrons. 
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Appendix 

In this appendix we discuss in some details the role of the evanescent operators in the 
matching procedure. Because of the IR divergences, the matrix elements of renormal- 
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ized evanescent operators do not vanish in the D ^ A hmit. Therefore, these operators 
contribute, at the NLO, to the matching of the physical operators. 

For the sake of simphcity, we consider here the abehan case. The inclusion of colour 
factors is trivial. 

Let us consider the calculation of the imaginary part of the double insertion of the 
operator 

Q = ibc)v-Aiud)v-A (56) 

in the left diagram of fig. |^. Neglecting the charm quark mass, the contribution to the 
amplitude in the full theory, proportional to the operators with the flavour structure bbdd, 
is 



7^(0) 



271 



1 e 1 , 

\ \ log At 

24 9 24 ^ml 



15 1 /i^ 
^ -7^^ + tt:^ log 



12 36 



12 



mt 



(57) 



where C(/i) is the AB = 1 Wilson coefficient of Q and 'J^n = 7^(1 ± 75)- Using the 
equations of motion, the result can be written as 
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[irMd^^Ld)- 
(b^^Sidl^Ld) + (5^ 
(67^7-7^6) (J7,7„7^^d) + 
{h'l''Yiib){d-f.7a7,Ld)] . 



We then introduce two physical operators and two evanescent operators (in the non-abelian 
case there are four of them, see eqs. ( |l4|) and (^0|) ) defined as 



01 = 67^6 d-ff^id , El = Iq^-ialuLh dYl^lid - 4 Oi , 

02 = b'j'^b d'j^Ld , E2 = hiilaluRb dYl"lLd - 16 O2 . 



(59) 



Using eqs. (^) and (|59D , the LO matching is easily performed and gives 



^(0) 
spect 



Gl\Vcb\^ml 



2n 



Cif^y E (co. Ok + CE, Ek 



k=1.2 



(60) 



''We explicitly checked that alternative definitions of evanescent operators, which differ from those in 
eq. (|9|) by terms of 0{e), give the same results in the final matching. 
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where the coefficient functions at the LO, including terms of 0{e), are 



1 13 1 /^^ 1 

'^^ = -3--8'~r''^^,^ ^-^ = -18' .g,. 
1 5 1 , /i^ 1 
6 18 6 mf 72 

The renormalized evanescent operators have to be inserted in the effective theory at 
the NLO. Due to the presence of IR poles, they acquire non-vanishing matrix elements in 
the D — >■ 4 limit. These finite contributions enter the matrix s in eq. (|17|) and, according 
to eq. ([T9|), contribute to the final determination of the Wilson coefficients at the NLO. 

A finite contribution is also obtained from the terms of 0{e) in the coefficient functions 
coj, of the physical operators. Indeed, in eq. (p!9D, these coefficients multiply the matrix s 
which, because of the IR divergences, contains poles in 1/e. Note that these terms also 
depend on the renormalization scale fi and are needed for reconstructing the proper UV 
scale dependence of the Wilson coefficients at the NLO. 
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